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Using Langevin dynamics simulations, we investigate tlie influence of polymer-pore interactions on 
the dynamics of biopolymer translocation through nanopores. We find that an attractive interaction 
can significantly change the translocation dynamics. This can be understood by examining the 
three components of the total translocation time r ~ n + r2 + rs corresponding to the initial filling 
of the pore, transfer of polymer from the cis side to the trans side, and emptying of the pore, 
respectively. We find that the dynamics for the last process of emptying of the pore changes from 
non-activated to activated in nature as the strength of the attractive interaction increases, and rs 
becomes the dominant contribution to the total translocation time for strong attraction. This leads 
to a new dependence of r as a function of driving force and chain length. Our results are in good 
agreement with recent experimental findings, and provide a possible explanation for the different 
scaling behavior observed in solid state nanopores vs. that for the natural a-hemolysin channel. 

PACS numbers: 87.15.A-, 87.15.H- 



I. INTRODUCTION 

The controlled transport of polymer molecules through 
a nanopore has received increasing attention due to its 
importance in biological systems and its potentially rev- 
olutionary technological ai^glications^ [iL [31. There is 
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studies devoted to this subject. In an important ex- 
periment, Kasianowicz et al. [l| demonstrated that an 
electric field can drive single-stranded DNA and RNA 
molecules through the water-filled a-hemolysin channel 
and that the passage of each molecule is signaled by 
a blockade in the channel current. These observations 
can be used to directly characterize the polymer length. 
Similar experiments have been done recently using solid 
state nanopores with more pre cisely controlled dimen- 
sions p, M, M m, m, sin. Hi, [13, (H- Currently, 
extensive effort is being taken to unravel the dependence 
of the translocation time r on the syst em p arameters such 
as the polymer chain len gth TV [ a, [ a , [2^ . [2a [2a . [27[ . [2^ 
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Meller et al. [1, [1] have shown how several differ- 
ent DNA polymers can be identified by a unique pat- 
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tern in an "event diagram". The event diagrams are 
plots of translocation duration versus blockade current 
for an ensemble of events. Patterns for a given poly- 
mer can be characterized uniquely by the blockade cur- 
rent, the translocation time and its distribution. Be- 
cause each type of polynucleotide gives rise to specific 
values of these three parameters, DNA molecules which 
differ from each other only by sequence can be distin- 
guished. At room temperature striking differences were 
found for the translocation time distributions of poly- 
deoxyadcnylic acid (poly(dA)ioo) and polydcoxycytidylic 
acid (poly(dC)ioo) DNA molecules. The translocation 
time of poly(dA) is found to be much longer, and its 
distribution is wider with a longer tail compared with 
the corresponding data for poly(dC). Moreover, the dif- 
ferences between the translocation behavior are accen- 
tuated at lower temperature. The origin of the differ- 
ent behavior was attributed to stronger attractive in- 
teraction of poly(dA) with the pore. Recently Krasil- 
nikov et al. [67| have investigated the dynamics of single 
neutral poly (ethylene glycol) (PEG) molecules in the 
a-hemolysin channel in the limit of a strong attractive 
polymer-pore attraction. The result for the residence 
time in the channel shows a novel non-monotonic behav- 
ior as a function of the molecular weight. The other ex- 
perimental data that point to the possible essential role of 
the monomer-pore interaction concerns the different con- 
flicting values of scaling exponents of r with N and with 
the applied voltage as reported in recent experiments. 
A linear dependence r iV was observed for polymer 
translocation through a-hemolysin channel P, [^ , while 
another experiment re por ted that r ~ 7V^-^^ w N'^'^ for 
a synthetic nanopore [2l|, where ly is the Flory expo- 
nent [H, [63|- As to the dependence of the transloca- 
tion time on the applied voltage for a-hemolysin channel, 
an inverse linear behavior [l[ is observed for polyuridylic 
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acid (poly(U)) while an inverse quadratic dependence [5[ 
is found for polydeoxyadenylic acid (poly(dA)). One pos- 
sible explanation for all these conflicting data is that the 
polymer-pore interaction depends crucially on the details 
of the pore structure (a-heniolysin channel vs synthetic 
nanopore) in addtion to being base pair specific. 

To date, most of the theoretical studies of the translo- 
cation of biopolymers through nanopre are based on mod- 
els in which the wall of the pore only plays a passive 
role in confining the polymer to the inside of the pore. 
There are only a few theoretical studies of such inter- 
action effects. Based on a Smoluchowski equation with 
a phenomenological microscopic potential to describe 
the polymer-pore interactions, Lubensky and Nelson [2^ 
captured the main ingredients of the translocation pro- 
cess. However, when comparing with experiments, their 
model is not sufficient. Numerically, Tian and Smith (60j 
found that attraction facilitates the translocation process 
by shortening the translocation time, which contradicts 
experimental findings [1, @- In a recent letter [s^, we 
used Langevin dynamics (LD) simulations to investigate 
the influence of polymer-pore interactions on transloca- 
tion. We found that with increasing attraction, the his- 
togram for the translocation time t shows a transition 
from Gaussian distribution to a long-tailed distribution 
corresponding to thermal activation over a free energy 
barrier. The N dependence of the entropic force leads to 
both the translocation time and the residence time in the 
pore showing a non-monotonic behavior as a function of 
N for short chains in the strong attraction limit. These 
results are in goo d agreement with the above experimen- 
tal data [l,[i,l63. 

In the present work, we further show that strong 
polymer-pore interactions can directly affect the effec- 
tive scaling exponents of t both with N and with the 
applied voltage, which provides a possible explanation 
for the different experimental findings [l|, |^, [2l| on these 
physical quantities. We provide a microscopic under- 
standing of how strong polymer-pore interaction influ- 
ences the translocation dynamics. This is done through 
analyzing the three quantities ti, T2 and T3 correspond- 
ing to initial filling of the pore, transfer of the polymer 
from the cis side to the trans side, and finally emptying 
of the pore, respectively. We find that the final process 
of emptying the pore T3 involves an activation barrier 
and completely dominates the translocation time in the 
strong attractive interaction limit. This leads to a strong 
dependence of the effective scaling exponents associated 
with the translocation time both on the strength of the 
attraction and the driving force. In addtion, we examine 
the waiting time and residence time distributions. These 
quantities are related to the translocation time but the 
waiting time provides more detailed information about 
the translocation dynamics, while the residence time is 
the more relevant quantity for direct comparison with 
with the experimental observation. 
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FIG. 1: A schematic representation of the system. The pore 
length L = 5 and the pore width = 3 (see text for the 
units) . 



II. MODEL AND METHODS 

In our numerical simulations, the polymer chains are 
modeled as bead-spring chains of Lennard-Jones (LJ) 
particles with the Finite Extension Nonlinear Elastic 
(FENE) potential. Excluded volume interaction between 
monomers is modeled by a short range repulsive LJ po- 
tential: ULjir) = 4e[(f - iff] -h e for r < and 
for r > Here, a is the diameter of a monomer, 

and e is the depth of the potential. The connectivity 
between neighboring monomers is modeled as a FENE 
spring with UpENEi'T) = —\kRQhi[l — v'^/Rq), where r 
is the distance between consecutive monomers, k is the 
spring constant and Rq is the maximum allowed separa- 
tion between connected monomers. 

We consider a 2D geometry as shown in Fig. [U 
where the wall in the y direction is described as sta- 
tionary particles within a distance a from each other. 
The pore of length L and width W in the center of 
the wall is composed of stationary black particles. Be- 
tween all monomer-wall particle pairs, there exist the 
same short range repulsive LJ interaction as described 
above. The pore-monomer interaction is modeled by a LJ 
potential with a cutoff of 2.5a and interaction strength 
Epm- This interaction can be either attractive or repul- 
sive depending on the position of the monomer from 
the pore particles. In the Langevin dynamics simula- 
tion, each monomer is subjected to conservative, fric- 
tional, and random forces, respectively, with [70| mf ; = 
-V([/lj + Ufene) + Fcxt - -I- Ff , where m is the 
monomer's mass, ^ is the friction coefficient, is the 
monomer's velocity, and is the random force which 
satisfies the fluctuation-dissipation theorem. The exter- 
nal force is expressed as Foxt = Fx, where F is the exter- 
nal force strength exerted on the monomers in the pore, 
and a; is a unit vector in the direction along the pore axis. 
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FIG. 2: Three components of the translocation process. 



In the present work, we use the LJ parameters e and 
a and the monomer mass m to fix the energy, length 
and mass scales respectively. Time scale is then given by 
tLj = (mcr^/e)^/^. The dimensionless parameters in our 
simulations are i?o = 2, fc = 7, ^ = 0.7 and kgT = 1.2 un- 
less otherwise stated. For the pore, we set L = 5 unless 
otherwise stated. A choice of = 3 ensures that the 
polymer encounters an attractive force inside the pore. 
The driving force F is set between 0.5 and 2.0, which 
correspond to the range of voltages used in the experi- 
ments [l|, Q . The Langevin equation is integrated in time 
by a method described by Ermak and Buckholz [7l[ in 
2D. Initially, the first monomer of the chain is placed in 
the entrance of the pore, while the remaining monomers 
are under thermal collisions described by the Langevin 
thermostat to obtain an equilibrium configuration. Typ- 
ically, we average our data over 2000 independent runs. 



III. RESULTS AND DISCUSSION 



A. Translocation time, waiting time and residence 

time 



The translocation time is defined as the time interval 
between the entrance of the first segment into the pore 
and the exit of the last segment. We can break down the 
translocation process into three components, as shown in 
Fig. 2. The total translocation time r can be written as 
a sum of three contributions r « ti +T2 + T3, where ri, T2 
and T3 correspond to initial filling of the pore, transfer of 
the polymer from the cis side to the trans side, and finally 
emptying of the pore, respectively. To shed light on the 
detailed translocation process, we examine the number of 
translocated monomers ritrans as a function of the time 
for a typical successful translocation event for N = 128, 
and two values of the monomer attractive interaction 




FIG. 3: Number of translocated monomers ntrans as a func- 
tion of the time for Spm = 1.0 and £^,71 = 3.0 under the 
driving force F = 0.5. For both strong and weak attraction 
strengths, n ^ r. For weak attraction strength Spm ~ 1.0, 
we find rg <^ T2 and thus t ~ T2. 



strength. The value £p„i = 1.0 corresponds to a weak in- 
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teraction whereas Spm = 3.0 corresponds to the strong at- 
traction limit. Here, ntrans = before the first monomer 
exits the pore and ntrans = N after the last monomer has 
threaded through the pore. As shown in Fig. 3, under 
the weak driving force _F = 0.5, ti is not sensitive to the 
attraction strength and ti <C T2 . T2 for the strong attrac- 
tion with Ejym ~ 3.0 is roughly twice as that for the weak 
attraction with Epm = 1.0. However, T3 depends strongly 
on the attraction strength. For Spm = 1.0, T3 <^ T2 and 
is basically negligible for the pore length L ~ 5. For the 
strong attraction limit with Spm = 3.0, the situation is 
totally different with T3 more than an order of magnitude 
larger than T2 , completely dominating the total contribu- 
tion to the translocation time. From Fig. 3, it can be 
seen that the number of translocated monomers oscillates 
around ntrans ~ 122, which corresponds to the beginning 
of the last stage of the translocation process, namely the 
emptying of the pore. This is due to the activated na- 
ture of the translocation process with a free energy dif- 
ference of AF = L{epm — F/2 — f{N)) between the final 
and the initial state. The term f{N) here accounts for 
the entropic driving force which should kick in at larger 
values of A^, and eventually saturate for very long poly- 
mers. This leads to the long oscillation time of the last 
few monomers with repeated forward and backward mo- 
tions. The final emptying of the pore corresponds to a 
rare crossing of the barrier. 

To provide more microscopic details of the transloca- 
tion process, we investigate the waiting time distribu- 
tion for different chain lengths N in the strong attrac- 
tion limit. The waiting time of monomer s is defined as 
the average time between the events that monomer s and 
monomer s+l exit the pore. In our previous work [46l.l47j 
for pure repulsive monomer-pore interactions, we found 
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FIG. 4: Waiting time of different chain lengths for Epm. = 3.0 
and F = 0.5. 



that the waiting time depends strongly on the monomer 
positions in the chain. For relatively short polymers, the 
monomers in the middle of the polymer need the longest 
time to exit the pore. Here, the waiting time of different 
chain lengths for Epm = 3.0 and F = 0.5 are shown in 
Fig. 4. ft can be seen that it takes much longer time for 
last three monomers to exit the pore, which is completely 
different from that for for pure repulsive monomer-pore 
interactions. This behavior correlates with the oscilla- 
tion of the last monomers as shown in Fig. 3. Here we 
should mention that due to the entropic factor f{N) in 
the barrier the waiting time for these last few monomers 
actually decreases in the range « 14 — 32 before satu- 
rating and even increasing slightly with further increase 
of iV. 

Under zero and low driving forces, the translocation 
probability is very small in the sense that many translo- 
cation events, once started do not finish all the way. In- 
stead, the polymer returns and exits to the cis side again. 
This means that the ti process of filling the pore do not 
get completed and the real translocation process corre- 
sponding to T2 and Ta never even get started. We define 
an additional residence time as the weighted average 
of the translocation time for the completed events and 
the return time for the events that start and return via 
the cis side. The significance of this quantity is that 
it corresponds to the experimentally measured average 
blockage time of the polymer in the nanopore which does 
not distinguish return events from the completed translo- 
cated events. For zero or low driving force {F < 0.5), 
the residence time is almost completely dominated by 
return events. We have calculated the residence time 
for Epm — 2.5 and 3 in Fig. 5. As shown in Ref. [sol . 
in the strong attraction case with Epm = 3.0, the N de- 
pendence of the residence time here is non-monotonic. 
This result of is in good agreement with experimental 
data of Krasilnikov et al. [13] where the residence time of 
a neutral PEG molecule in a-hemolysin pore was mea- 
sured. Here, we further show that for Epm = 2.5, 
increases with increasing N . It indicates that the strong 
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FIG. 5; Residence time Tt as a function of the chain length 
for Epm = 3.0 and Epm ~ 2.5 under the driving force F — 0. 
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attraction plays an essential role in the observed non- 
monotonic behavior. 



For 



3.0, the distribution of is shown in Fig. 



6. One obvious feature is the existence of two groups. 
The first group with shorter corresponds to the events 
where one end of the chain accesses the pore, and then 
quickly returns back. For the second group with longer 
Tr, the residence time is still about 99.8% due to return 
events for Epm = 3.0. In the strong attraction limit, once 
the attractive force reaches its maximum when the pore 
is fully filled by monomers, it takes a very long time for 
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FIG. 7: Translocation time as a function of the temperature 
for both strong and weak attraction strengths (e^m = 3.0 and 
1, respectively) under the driving force F = 0.5. The chain 
length N = 128. 
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FIG. 8: Translocation time as a function of the driving forces 
for both strong and weak attraction strengths, Epm = 3 and 
1. The chain length iV = 128. 



the polymer to return back due to frequent backward and 
forward events. 



B. Dependence of translocation time on various 
parameters 

1. Translocation time as a function of temperature 

Fig. 7 shows the translocation time r as a function of 
the temperature for different attraction strengths. For 
the whole examined range of temperatures, r decreases 
very slightly with increasing temperature for a weak at- 
tractive strength of Epm = 1-0. However, for the strong 
attractive strength Spm = 3.0, with increasing tempera- 
ture T first rapidly decreases and then approaches satu- 
ration at higher temperatures. At higher temperatures, 
the differences between translocation times for weak and 
strong attractive strengths become very small. This tem- 
perature dependence of translocation time is in good 
agreement with experiments 



2. Translocation time as a function of the driving force 

In the weak attraction (i.e. non-activated) region, the 
overall r is determined mainly by T2 and its dependence 
on the driving force scales as F~^. This simple scaling 
behavior can be understood by considering the steady 
state of motion of the polymer through the nanopore. 
The average veleity is determined by balancing the fric- 
tional damping force (proprtional to the velocity) with 
the external driving force. This leads to an average ve- 
locity proportinal to the driving force F, and hence a 
translocation time r ~ F^^. In Fig. 8 we show the 
dependence of the translocation time r on the driving 
force. It can be seen that in the weak interaction limit 



for Spm = 1.0 the data are very close to the linear scaling 
behavior r ~ F^^ as predicted. For strong attractive 
intcaction with Epm — 4.0, the situation is more compli- 
cated. For weak driving forces (F < 2), one is in the the 
activated region where the inverse of the translocation 
time obeys an Arrhenius form. However, the driving force 
F affects both the activation barrier and the prefactor, 
leading to a complicated dependence of t on the driving 
force that does not have a simple power law scaling form 
as seen in Fig. 8 for the Epm = 3.0 result. Insistence 
on fitting the data with a power law scaling form will 
lead to an effective scaling exponent that changes with 
the value of the driving force. Finally, beyond a criti- 
cal force, the activation barrier disppears and one should 
obtain asymptotically the t ~ F~^ behavior just as in 
the weak interaction case. This whole scenario is very 
similar to the sliding friction of an adsorbed layer under 
an external driving force [73 | . 

The above theoretical considerations lead to a possible 
explanation of recent apparently conflicting experimental 
data. Polyuridylic acid (poly(U)) has a weak interaction 
with the pore, and it is not surprising that an inverse 
linear dependence of the translocation time on applied 
voltage was observed in experiments on the translocation 
of poly(U) However, Poly(dA) has much stronger in- 
teraction with the pore compared with poly(U). Thus 
it should be in the strong interaction activated region 
with a larger effective scaling exponent. Indeed, an in- 
verse quadratic dependence of the translocation time on 
applied voltage had been experimentally observed for 
poly(dA) 0. In view of our theoretical considerations, 
it would be desirable to have measurement over a larger 
range of the applied voltage to see the predicted change 
of effective scaling exponent. 
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FIG. 9: Translocation time as a function of the chain length 
for Epm = 3.0 under F — 1.0 and F = 2.0, respectively. 

3. Translocation time as a function of chain length 

Previously, wc have established that for pure repul- 
sive polymer-pore interactions, the dependence of the 
translocation time on the length of the polymer scales 
as r ~ N^'^ for N < 200 and crosses over to a new scal- 
ing regime r ~ N^^" for larger values of N pEli?!. [iot. 
In the presence of weak interaction between the monomer 
and the pore, the qualitative dependence on the length 
of the polymer remains the same. For stronger attrac- 
tive strength Spm — 3.0, the scaling exponent of r with 

for 64 < iV < 400 becomes strongly dependent on the 
driving force, with no indication of crossover behavior as 
shown in Fig. 9. We find r ~ N^-^^ for F = 2.0, which 
is close to T ^ N'^'^ with the Flory exponent i' = 0.75 in 
2D [el, m, and t - N^-^'' for F = 1.0. The novel de- 
pendence on the length of polymer is due to the change 
from the non-activated regime for weak attractive or pure 
repulsive interaction to an activated regime for strong at- 
tractive interaction. 

Experimentally, a linear dependence t ^ N was ob- 
served in experiments [l|, Q for polymer translocation 
through a-hemolysin channel, in contrast to the r ~ N'^'^ 
scaling observed for pol ymer translocation through the 
solid-state nanopore |2l| . This difference can be under- 
stood in light of our present results concerning the in- 
fluence of the different polymer-pore interaction on the 
length dependence of the translocation time. For a syn- 
thetic pore, there is at most a very weak attractive inter- 
action between the polymer and the pore and one expects 
the scaling behavior r ~ N'^" to hold for N < 200. How- 
ever, a stronger attractive interaction is expected to exist 
between the different bases and the a-hemolysin channel. 
For the models studied in this work, it changes the scal- 
ing behavior from r ~ N^'^ to t ^ N. This provides a 
possible explanation for the difference of the experimen- 
tal observations in the different nanopores P, H, [2l[ . 




N 

FIG. 10: Translocation time r as a function of the chain 
length for Epm = 3.0 and Epm = 2.5 under the driving force 
F = 0.5. 



Under a strong attractive force with Epm = 3.0 and 
a weak driving force F = 0.5, the translocation time r 
has a qualitatively different dependence on N as com- 
pared with the pure repulsive or weak attractive pore 
interaction. Here we should mention that for F = 0.5 we 
cannot access N > 128 as the translocation time becomes 
too long to be feasible for numerical comuputation. As 
shown earlier in Ref. [13 and here in Fig. 10, the translo- 
cation time displays novel non-monotonic behavior with 
a rapid increase to a maximum at A'^ ~ 14, followed by 
a decrease for 14 < < 32 and an increases again for 

> 32. The eventual increase in the large A^ limit is due 
to the T2 contribution for longer chains. The observed 
non-monotonic behavior is also reflected qualitatively in 
the waiting time distribution as shown in Fig. 4. As 
shown in Fig. 10, with decreasing an attractive force to 
Epm = 2.5, this non-monotonic behavior vanishes. 

To understand the microscopic origin, in Fig. 11 we 
show Ti -f- r2 as a function of the chain length for differ- 
ent attraction strengths under the driving force F = 0.5. 
For 32 < Af < 200, Ti T2 ~ N^" is observed, irrespec- 
tive of attraction strengths. This indicates that the novel 
non-monotonic behavior shown in Fig. 10 in the strong in- 
teraction limit is again due to the pore-emptying process 
corresponding to dominating the translocation time in 
strong interaction limit. 
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FIG. 11: Ti +r2 as a function of the chain length for different 
Epm under the driving force F = 0.5. 

IV. CONCLUSIONS 

111 this work, we have studied the dependence of 
the translocation time on the temperature, attraction 
strength, driving force and the chain length. To ana- 
lyze the influence of the attractive interaction in more 
detail, we have considered the three components of the 



translocation time t m ti + T2 + t^, which were exam- 
ined as a function of the attraction strength. Here Ti, 
T2 and Ta correspond to initial fllling of the pore, trans- 
fer of polymer from the cis side to the trans side, and 
emptying of the pore, respectively. We find that ti <C T2 
for both weak and strong attraction strengths, for TV in 
the typical range used in the experiments. However, 
is sensitive to the presence of an attractive interaction 
and changes from a value much less than T2 for weak 
interactions to the dominant contribution to the overall 
translocation time due to the rare activated event nature 
of the final emptying of the pore. This leads to a drastic 
change of the translocation dynamics and various scaling 
exponents as a function of the strength of the attractive 
monomer pore interactions. Our theoretical results are in 
good agreement with recent experimental data 0, 0, [67j . 
They also provide a possible explanation for the differ- 
ence of the scaling behaviors with regard to the driving 
force and the length of polymers observed using different 
types of nanoporesfll. [sl I2H . 
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